Aggregation in a mixture of Brownian and ballistic wandering particles 
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In this paper, we analyze the scahng properties of a model that has as limiting cases the diffusion- 
limited aggregation (DLA) and the ballistic aggregation (BA) models. This model allows us to 
control the radial and angular scaling of the patterns, as well as, their gap distributions. The par- 
ticles added to the cluster can follow either ballistic trajectories, with probability Pba, or random 
ones, with probability Prw = 1 — Pba- The patterns were characterized through several quantities, 
including those related to the radial and angular scaling. The fractal dimension as a function of Pha 
continuously increases from d/ ~ 1.72 (DLA dimensionality) for Pba = to d/ ~ 2 (BA dimension- 
ality) for Pba = 1. However, the lacunarity and the active zone width exhibt a distinct behavior: 
they are convex functions of Pba with a maximum at Pba ~ 1/2. Through the analysis of the an- 
gular correlation function, we found that the difference between the radial and angular exponents 
decreases continuously with increasing Pba and rapidly vanishes for Pba > 1/2, in agreement with 
recent results concerning the asymptotic scaling of DLA clusters. 

PACS numbers: 61.43.Hv,05.40.Fb,47.53.+n,47.54.-r 



I. INTRODUCTION 



Nonequilibrium growth models of clusters of identical 
particles have attracted a lot of attention in the last two 
decades [H 0- Certainly, the most studied one is the 
diffusion-limited aggregation (DLA) model proposed by 
Witten and Sander [3] for the modeling of metal par- 
ticle aggregation. In spite of its simplicity, this model 
is a noteworthy example in which a very simple algo- 
rithm generates disordered patterns with nontrivial scal- 
ing 0, 0, S H- I^ the DLA model, the particles are 
released at points distant from the cluster and follow ran- 
dom walks of unitary steps, which are stopped when these 
particles contact the cluster and irreversibly aggregate to 
it. Another important aggregation model with nontrivial 
scaling is the ballistic aggregation (BA) 0, 13 , in which 
ballistic trajectories at random directions replace the ran- 
dom walks characteristic of the DLA model. Differently 
from the DLA case, the patterns generated by the BA 
model are asymptotically homogeneous (fractal dimen- 
sion equal to the space dimension), but this asyrnptotic 
regime is characterized by a power law approach j^ [l^l . 

Several generalized versions of the DLA model were 
proposed |l|, including those concerned with the tran- 
sition between DLA and BA models [13, lU H^ ■ In ^ 
recent work, Ferreira et al. [lOJ built a model in which 
the particles follow random trajectories with drifts at ran- 
dom directions. The bias is controlled by a parameter A 
that, when varied, leads the model from BA (A = 0) to 
DLA (A=l). The clusters are asymptotically non- fractal 
(exhibiting the BA-like scaling with fractal dimension 
df = 2) for any A 7^ 1, but following a DLA- like scaling 
{df K, 1.71) for short length scales. Using the concept of 



scaling functions, these ideas can be summarized as 

M{1) = f^-^JillO, (1) 

in which 



/(^) 



const, if a; <IC 1 

^d-doLA if x > 1. 



(2) 



In these equations, I is the size of the cluster, d is the 
space dimension, dj^^A is the DLA fractal dimension, 
and ^ is the characteristic crossover radius from DLA- 
to BA-like scaling regimes. The last diverges as ^ ~ 
1 — \\~'^ where v — 0.61 ±0.01. This transition between 
DLA and BA regimes was introduced by Meakin [13 in 
a early generalization of the DLA model. This model 
is very similar to that analyzed by Ferreira et al. |10J . 
The main difference is that the drift of all trajectories 
is fixed at a lattice direction in the Meakin's model, in 
which, along the trajectory, the particle is moved one 
step in the drift direction with probability P, or moves 
at random to one of its four nearest-neighbor sites with 
probability 1 — P. The model generates patterns with 
a growth tendency along the opposite direction of the 
drift. In this work, Meakin argued that the crossover 
radius scales as ^ = P~^ while Nagatani 14], through 
a real-space renormalization group approach, found ^ ^ 

p-i/{d-df)_ 

Early studies indicated that the DLA aggregates ex- 
hibit distinct scaling exponents for radial and angular 
correlations. The two-point correlation function F, which 
gives the probability to find a particle at a distance be- 
tween r and r -f 5r, Sr <C r, from another particle, is 
given by [l| 

F(r)~r-«^ (3) 



where a^ = d ~ df. For ofF-lattice DLA, the exponent 
found was a^ ~ 0.29. In turn, the angular correlation 
function $ij, which gives the probability to find two par- 
ticles with angular separation 9 measured in relation to 
the origin at a distance R, follows a power law decay for 
small 9 



<^r{0) 



(4) 



For uniform self-similar fractals a,. = ag, but this was 
not observed for the DLA model, for which ag ~ 0.41 

m 

In this work, we used careful simulations to analyze 
the transition between DLA- and BA-like regimes in a 
model where both random and ballistic trajectories are 
included. Differently from the previously described mod- 
els [i0,[i3; the trajectories are not characterized by dis- 
tinct behaviors at short (random) and large (ballistic) 
scales. Indeed, we are interested in the scaling of clusters 
grown from a mixture of particles, which can follow ci- 
ther ballistic or random trajectories. The paper outline 
is the following. In section^ the algorithm and the com- 
puter implementations are described. In section UllI the 
quantities describing the scaling of the cluster are intro- 
duced and the corresponding results discussed. Finally, 
we summarize the results of the paper in section Hvl 



II. MODEL AND COMPUTER 
IMPLEMENTATION 

The present model is concerned with the aggregation 
of particles which can follow either ballistic or random 
trajectories. Notice that there are two kinds of inde- 
pendent trajectories, which do not exhibit both random 
and ballistic components simultaneously, as those used 
in reference ,lQj . Two-dimensional off-lattice simulations 
starting from a single particle of unitary radius stuck to 
the origin were considered. The model's rules are very 
simple. For each time step, a kind of trajectory, Brown- 
ian or ballistic, is chosen with probabilities Prw and Pha, 
respectively. Then, particles are launched one at a time 
using the selected rule until one of them sticks to the 
cluster. Notice that in these rules we are taking in to 
account only particles, which are added to the cluster 
and not those that get too far from the cluster and are 
excluded. 

The particles are released at a launching circle of ra- 



dius ri > Vmax centered at the origin, where 



is the 



largest distance from the origin among all the aggregated 
particles. As in the off-lattice DLA model |13] , the cluster 
grows when a walker contacts any particle of the aggre- 
gate and sticks to the corresponding position. Also, if the 
walker crosses a killing circle of radius rt > n centered 
at the origin, it is discarded and a second one released at 
a new random point on the launching circle. The previ- 
ously introduced variables (ri and rk) should be as large 
as possible. However, computational limitations restrict 
their values. For sake of convenience, we introduce a new 
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FIG. 1: Schematic representation of the "optimized random 
trajectories". We show a DLA aggregate and a mesh of cells 
25m X 25i„. Long steps are forbidden in the gray boxes and 
allowed in the white ones. Also, two long steps are illustrated. 



variable E by the definition r; = r^aax + ^- The choices of 
b and r^ depend on the kind of trajectory. For random 
walks, Ti can be a few particle diameters larger than r^ax 
whereas r^ must be very large [l| . Consequently, we used 
(5 = 5 and r]~ = lOOr^aa:- For ballistic trajectories unde- 
sirable shadow effects are reduced as r; increases [l9H2(lj . 
Thus, we used b = lOOrmaa: + 800. Finally, ballistic tra- 
jectories which crosses the launching circle can not turn 
back inside the circle and, so, we chose r^ = r; -|- 2. 

Large scale off-lattice simulations with rigorous statis- 
tical sampling require very efficient algorithms. Firstly, 
it is necessary to improve the efficiency of the routines for 
the trajectories. In the case of random walks, we used a 
standard method, in which particles outside the launch- 
ing circle take long steps ^out if these steps do not bring 
up a particle inside the launching circle. An adequate 
choice is ^o^t = max(r — Vmax — 5, 1), where r is the dis- 
tance of the walker from the origin. Also, in the inner 
region of the circle delimiting the cluster there are large 
empty regions. So, we adopted an algorithm like those 
used by Ball and Brady f2l| , in which particles inside the 
launching circle can take a long step of length ^i„ if they 
cannot cross any part of the aggregate along this long 
step. 

An efficient determination of ^m is decisive in order 
to guarantee the success of the procedure. In order to 
accomplish this, we define a coarse-grained mesh with 
cells of size 2fi„ x 2^i„ as illustrated in figure ^ The 
boxes depicted in gray are those in which the random 
walk could reach the cluster after a step of length ^i„ 
and, consequently, a long step is forbidden. In this case, 
we have two options: the particle executes a unitary step 



or tries a shorter step of length ^^'^ , where 1 < S,ln < iin , 
using another auxiUary coarse-grained mesh. Indeed, we 
can use several auxiliary meshes in order to maximize 
the algorithm efficiency. In this work, three values, ^i„ — 
8, 16 and 32, were used. Through this procedure, large 
simulations become up to two orders of magnitude faster 
than those using only long steps outside the launching 
circle. 

Finally, an efficient search mechanism for determining 
when the walker has or not a contact with the aggregate 
is indispensable. We used two procedures for this task. 
In the first one, the particle positions are mapped on a 
square lattice by the approximation of their real coordi- 
nates to the nearest integer, producing an on-lattice aux- 
iliary cluster. So, the verification of the contact is done 
only if the particle is inside a nearest or a next-nearest 
empty site neighboring the on-lattice auxiliary cluster. 
In the second optimization procedure, a coarse-grained 
mesh of cells 8x8 was used to limit the verification on 
a region near the walker position. 



III. 



RADIAL AND ANGULAR SCALING 



The model exhibits a smooth transition from DLA to 
BA-like growth patterns as the parameter P^a varies from 
to 1 . In figure [3 in which growth patterns for dis- 
tinct probabilities Pba are shown, one can observe clusters 
that are qualitatively very similar to the DLA aggregates 
for Pba < 0.5, whereas BA-like clusters are observed for 
Pba > 0.95. These patterns were grown until a part of 
them reaches the edge of a region of size L x L, where 
L = 750. Notice that Pta = 0.5 is an important case, in 
which a half of the particles were grown using DLA model 
and the another half were grown using the BA model. 
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FIG. 2: Growth patterns for distinct values of the parameter 
Pba-- (a) 0.1, (b) 0.3, (c) 0.5, (d) 0.7, (e) 0.9, (f) 0.95, and (h) 
0.99. All simulations were stopped when the cluster reaches 
the border of a region 750 x 750 centered at the origin. The 
number of particles varies from A'^ ~ 3 x 10* (Pta ~ 0.1) to 
iV ~ 1.6 X 10^ {Pba = 0.99). 

The scaling properties of the clusters were character- 
ized by the radial and angular distributions of the parti- 
cles, as well as, the size of their active (growing) zone. We 
grew 2d off-clusters limited by a region of size 10^ x 10^. 



Again, the simulations were stopped when the cluster 
reaches border of this region. The radial scaling was 
characterized by the fractal dimension and the lacunar- 
ity 1]. The fractal dimension df and the lacunarity A 
were determined using the mass-radius and the radius of 
gyration methods. The mass-radius method consists of 
counting the number of particles M in a hypersphere of 
radius r centered at the origin. We expect a power law 
relation for self-similar clusters 



Af(r) 






(5) 



where the amplitude Am is a measure of the lacunarity. 
To be more precise, the DLA clusters are characterized by 
a multiscaling, which leads to a slow convergence to the 
asymptotic scaling exponents. However, large scale sim- 
ulations confirm that these clusters are asymptotically 
fractal. For details see references [l|,|22. The radius of 
gyration Tg is defined as 



N^ 



(6) 



where r^ is the distance from the origin of the zth aggre- 
gated particle and N the total number of particles in the 
cluster. The gyration radius scales with the number of 
particles as 



Tn = BiV^ 



(7) 



where v = 1/df and B can be related to the lacunarity 
trough the relation 



A J- 



(8) 



Notice that Am ^ Ag since their definitions are different. 
But, they exhibit the same qualitative behavior because 

Am OC Ag. 

In figure |3Ja), the double logarithm plots of the curves 
M{r) against r are shown for several Pba values. The 
slope continuously increases from df = 1.723 ± 0.007 for 
Pba = 0.10 (a value close to the well-known DLA frac- 
tal dimension df = 1.715 ± 0.004 ^) towards 2.0 as 
Pba increases, as can be seen in figure |2Ib), in which the 
fractal dimension is plotted as a function of Pba ■ Remark- 
ably, instead of the qualitative behavior observed for the 
fractal dimension (an increasing function of Pba), the la- 
cunarity is a convex function of Pba with a maximum at 
Pba ~ 1/2, as shown in inset of figure |3Ia). These re- 
sults show that if the number of DLA particles is larger 
than the number of BA ones, an increasing Pba produces 
patterns with more dense branches (larger df values) but 
with larger empty regions (larger lacunarity). However, 
when the number of BA particles exceeds the number of 
DLA ones, the width of the branches increases while the 
size of empty regions decreases with Pba ■ An explanation 
for such a behavior can be obtained from the active zone 
analysis described next. 
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FIG. 3; (a) Mass against radius for distinct Pta values: 0.10, 0.20, 0.30, 0.40, 0.50 0.70, 0.90, and 0.99 from the bottom to 
the top. (b) Fractal dimension as a function of Pta determined from the mass-radius (•) and gyration radius (D) methods. 
In inset, we show the lacunarity as a function of Pta obtained from mass-radius (•) and radius of gyration (□) methods. For 
these simulations we used clusters of size L = 10^ and averages done over 10^ independent samples. The number of particles 
in these clusters varies from 7 x 10® for Pta = 0.1 to 3 x 10^ for Pta ~ 0.99. 



We define the active zone as the fraction / of the latest 
particles added to the cluster, and its width ct/ as the 
standard deviation of the distances from the origin of 
these particles. So, 



-? 



fN 



N 

E 

t=Nf 



where iV/ = (1 - f)N and 



(r) 



1 



in - {r)r 



N 

E 

l=Nf 



(9) 



(10) 



In all simulations we used / = 0.1. In figure 01 the 
width of active zone as function of the number of particles 
is shown for two distinct Pi,a values. For P^a < 0.5, 
the simulations indicate that there is a single power law 
regime a ~ N'^ while a crossover between two power 
laws a ~ N'^ for small length scales and a ^ N'' for 
large ones is observed for Pta > 0.5. In the asymptotic 
regime (TV — > oo), the exponent ^ — > i^ for p < 1 as 
observed in the DLA model [Tg- For Pf,a > 0.5, we 
found C' ~ 0.46, 0.42 and 0.36 for Pba = 0.5, 0.7 and 
0.9, respectively. Indeed, our simulation suggests that 
C -^ 0.33 as Pba -^ 1, in agreement with the statement 
that interface scaling of the BA belongs to the Kardar- 
Parisi-Zhang universality class [23 • To be more precise, 
the interface scaling of ballistic deposition models were 
widely studied and rigorous simulations confirmed that 
(T-7V1/3 0,0. 

A very interesting feature emerges from the behavior 
of cr as a function of P^a for a fixed mean radius of the 




FIG. 4: Active zone width cr as a function of the number of 
particles for Pta ~ 0.20 (•) and 0.90 (o). The straight lines 
correspond to power law fits. In inset, we show cr as a function 
of Pta for a mean active zone radius r ~ 340. The line is a 
cubic fit used as a guide to the eyes. For these simulations we 
used clusters of size L — 10^ and averages were done over 10'^ 
independent samples. 



active zone, as shown in inset of figure ^ The active 
zone width is a convex function of Pba with a maximum 
at Pba = 1/2, which explains the lacunarity dependence 
on Pba- The onset of this qualitative behavior can be 



understood through an analysis of the screening proper- 
ties of the clusters: for Pi,a ^ 0, the patterns have large 
grooves and the ballistic trajectories can reach the inner 
regions of the cluster forbidden for random walks. Thus, 
the screening is reduced in relation to the original DLA 
model and, consequently, the active zone is enlarged. As 
Pba increases, more particles break the screening enlarg- 
ing the active zone. However, if too many particles arrive 
at the inner regions of the clusters, the size of the their 
grooves and, consequently, their active zones are reduced. 
Therefore, the convexity of a versus Pba was elucidated. 
Intriguingly, this change of behavior seems to occur when 
the number of DLA and BA particles are equal, but we 
cannot provide a rigorous justification besides the numer- 
ical evidence. 




FIG. 5: Two point angular correlation functions evaluated at 
R = 300 for several values of Pta'. 0.1, 0.2, 0.3, 0.5, 0.7, and 
0.9 from the bottom to the top, respectively. The dashed lines 
correspond to the power law fits for small 0. In inset, the an- 
gular and radial correlation function exponents are compared. 
The size of the clusters is 5000 x 5000 and ^r was evaluated 
using R — 300 and SR = 15. The averages were done over 
N = 10^ independent samples. 



One of the most complex features of DLA clusters is 
their multiscaling structures |a, iWL Il7| . which was re- 
cently analyzed by Mandelbrot et al. at the light of the 
angular gap distributions of the clusters Q • They found 
that the angular and radial distributions of mass scale 
with the same exponent D « 1.71 only for very small an- 
gular gaps. This happens for asymptotically large clus- 
ters [M K, 10®) when the number of main branches in- 
creases leading to a more uniform space fulfillment. From 
figure [3 one can see that angular gap distributions vary 
with the parameter P^a- Therefore, a angular scaling 
analysis becomes necessary in order to obtain a more 
elaborate characterization of model. For this task, we 
followed the work of Meakin and Vicsek [I3, in which 



they analyzed the angular correlations of the DLA model 
using the following procedure. A narrow circular shell of 
radius R and width SR <^ R centered at origin is divided 
in K identical sectors of area RS95R, where S9 — 2tt/K. 
For a box centered at the polar coordinates {R,(p), the 
quantity pr{4>) = k is defined, where k is the number of 
particles within this box. The two-point angular correla- 
tion function is given by 



K/2 



$fl,(^) = ^J2p"(^ + nS0)pR{nSe). 



(11) 



The exponent obtained for the relatively small DLA clus- 
ters (~ 10^) using this procedure was ae « 0.41, which 
is definitely different from a,, w 0.29. 

Figure illustrates the angular correlation functions 
evaluated at _R = 300 for distinct Pba values. In all 
curves, ^r{0) decays as a power law for small angles fol- 
lowed by a global minimum and reaches a constant value 
for large angles. This behavior qual itatively reproduces 
the finds by Meakin and Vicsek [131, but the power law 
exponents decay with Pba as illustrated in inset of figure 
13 in which ae and a^, as well as, the difference between 
them Aa are plotted as functions of Pba- The expo- 
nent ar was determined using the relation a^ — 2 — df, 
where df is the mean fractal dimension obtained through 
the mass-radius and radius of gyration methods (figure 
Elb)). The difference Aa « 0.12 observed for the DLA 
model decreases with Pba- The minima of the curves 
are interesting quantities since they represent estimates 
of the half angular separation between two consecutive 
branches, the gaps studied by Mandelbrot et al. Q- The 
minima position shift to small angles and become less 
evident with increasing Pba- Indeed, the extrapolation 
of minimum positions suggests that their values vanishes 
at Pba = 1 following a power law 9,nin ~ (1 - Pba)"^, 
where •& ~ 0.45, as shown in figure The minima were 
determined using a cubic fit around their positions, as 
indicated in insets of figure El This power law relation 
means that the angular gaps are present, even that small, 
for any Pba 7^ 1. However, one can resolve from inset of 
figurcEl the difference Aa vanishing for Pba ~ 0.5. Thid 
behavior of Aa agrees with the concluding claim of Man- 
delbrot et al- [5J concerning the asymptotic scaling of or- 
dinary off-lattice DLA clusters. Indeed, in our model the 
gaps are reduced with increasing Pba and, consequently, 
we observed the same result of reference [5(| , but without 
considering simulations with more than 10® particles. 



IV. SUMMARY 

In this work, we analyzed an aggregation model in 
which two kinds of particles are considered: a given frac- 
tion following Brownian motions and the remaining fol- 
lowing ballistic trajectories. The model, which includes 
diffusion-limited aggregation (DLA) and ballistic aggre- 
gation (BA) models as limit cases, is controlled by a prob- 
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FIG. 6: Behavior of the minima presented in angular cor- 
relation functions for distinct probabilities Pf,a. The dashed 
Une is a power law fit corresponding to an exponent i9 ~ 0.45. 
Data in the range Pi,a £ [0.1, 0.95] are shown. The insets show 
the linear plots of ^r for the regions near to the minima for 
Pba = 0.7 (bottom) and Pba ~ 0.9 (top). The solid lines are 
the cubic fits used to determine the positions of the minima. 



ability Pba that correspond to the fraction of B A particles 
which constitutes the aggregate. Large scale simulations 
and careful statistical sampling were used to characterize 
the radial and angular scaling of off-lattice aggregates. 
Specifically, we measured the fractal dimension df, the 
lacunarity, the width of active (growing) zone cr, and the 
angular correlation function <!>/;• of clusters grown for dis- 
tinct Pba values. 

The fractal dimension df continuously increases from 



df = 1.715 ± 0.004 for Pba = (DLA fractal dimen- 
sion) to d/ « 2 (BA fractal dimension). However, the 
lacunarity and the width of the active zone for a fixed 
radius are convex functions of the probability Pba with 
maxima at Pba ~ 1/2. The convexity of the lacunarity is 
directly related to the convexity of the active zone width, 
which in turn, was explained by the screening impairing 
the wandering particles to reach the inner regions of the 
clusters. 

The angular correlation function $i?(6') evaluated at 
R = 300 exhibits a power law decay for small angles with 
a exponent that is a decreasing function of Pba- The dif- 
ference between angular and radial exponents observed 
for the original off-lattice DLA model [13 was found, but 
it rapidly vanishes for Pba > 1/2. Also, the curves ^ji 
versus 6 exhibit global minima representing the half an- 
gular separations between two consecutive branches. The 
simulations show that the position of the minima van- 
ishes as Pba -^ 1 following a power law 9min ~ (1 — Pba)^ 
with i) « 0.45. These results indicate that the patterns 
are asymptotically homogeneous when Pba — *■ 1. In sum- 
mary, in this work we show how the scaling properties 
of the aggregates can be controled by varying the con- 
centration of the kinds of wandering particles involved in 
the growth process. Finally, it is important to empha- 
size that these results are in agreement with the angular 
^p distribution analysis performed by Mandelbrot et al. 

, in which they found that the original DLA clusters 
exhibit the same asymptotic angular and radial scaling. 
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